    fileName path
    (
        UMean.rootPath()
       /UMean.caseName()
       /functionObjects::writeFile::outputPrefix
       /"graphs_My"
       /UMean.instance()
    );
    mkDir(path);

    scalarField UMeanXvalues
    (
        channelIndexing.collapse(UMean.component(vector::X)())
    );

    scalarField UMeanYvalues
    (
        channelIndexing.collapse(UMean.component(vector::Y)())
    );

    scalarField UMeanZvalues
    (
        channelIndexing.collapse(UMean.component(vector::Z)())
    );

    scalarField RxxValues(channelIndexing.collapse(Rxx));
    scalarField RyyValues(channelIndexing.collapse(Ryy));
    scalarField RzzValues(channelIndexing.collapse(Rzz));
    scalarField RxyValues(channelIndexing.collapse(Rxy, true));
    
    scalarField pPrime2MeanValues(channelIndexing.collapse(pPrime2Mean));
    
    scalarField T_0_001MeanValues(channelIndexing.collapse(T_0_001Mean));
    scalarField T_0_005MeanValues(channelIndexing.collapse(T_0_005Mean));
    scalarField T_0_025MeanValues(channelIndexing.collapse(T_0_025Mean));
    scalarField T_0_05MeanValues(channelIndexing.collapse(T_0_05Mean));
    scalarField T_0_125MeanValues(channelIndexing.collapse(T_0_125Mean));
    scalarField T_0_25MeanValues(channelIndexing.collapse(T_0_25Mean));
    scalarField T_0_71MeanValues(channelIndexing.collapse(T_0_71Mean));
    scalarField T_1MeanValues(channelIndexing.collapse(T_1Mean));
	
//add T no time mean
	scalarField T_0_001Values(channelIndexing.collapse(T_0_001));
    scalarField T_0_005Values(channelIndexing.collapse(T_0_005));
    scalarField T_0_025Values(channelIndexing.collapse(T_0_025));
    scalarField T_0_05Values(channelIndexing.collapse(T_0_05));
    scalarField T_0_125Values(channelIndexing.collapse(T_0_125));
    scalarField T_0_25Values(channelIndexing.collapse(T_0_25));
    scalarField T_0_71Values(channelIndexing.collapse(T_0_71));
    scalarField T_1Values(channelIndexing.collapse(T_1));
	
	
    
    scalarField T_0_001Prime2MeanValues(channelIndexing.collapse(T_0_001Prime2Mean));
    scalarField T_0_005Prime2MeanValues(channelIndexing.collapse(T_0_005Prime2Mean));
    scalarField T_0_025Prime2MeanValues(channelIndexing.collapse(T_0_025Prime2Mean));
    scalarField T_0_05Prime2MeanValues(channelIndexing.collapse(T_0_05Prime2Mean));
    scalarField T_0_125Prime2MeanValues(channelIndexing.collapse(T_0_125Prime2Mean));
    scalarField T_0_25Prime2MeanValues(channelIndexing.collapse(T_0_25Prime2Mean));
    scalarField T_0_71Prime2MeanValues(channelIndexing.collapse(T_0_71Prime2Mean));
    scalarField T_1Prime2MeanValues(channelIndexing.collapse(T_1Prime2Mean));
    
    scalarField tke_CkMeanValues(channelIndexing.collapse(tke_CkMean));
    scalarField tke_DkMeanValues(channelIndexing.collapse(tke_DkMean));
    scalarField tke_EpikMeanValues(channelIndexing.collapse(tke_EpikMean));
    scalarField tke_TkMeanValues(channelIndexing.collapse(tke_TkMean));
    scalarField tke_PkMeanValues(channelIndexing.collapse(tke_PkMean));
    scalarField tke_PikMeanValues(channelIndexing.collapse(tke_PikMean));
    
    scalarField uu_CkMeanValues(channelIndexing.collapse(uu_CkMean));
    scalarField uu_DkMeanValues(channelIndexing.collapse(uu_DkMean));
    scalarField uu_EpikMeanValues(channelIndexing.collapse(uu_EpikMean));
    scalarField uu_TkMeanValues(channelIndexing.collapse(uu_TkMean));
    scalarField uu_PkMeanValues(channelIndexing.collapse(uu_PkMean));
    scalarField uu_PikMeanValues(channelIndexing.collapse(uu_PikMean));
    scalarField uu_PhikMeanValues(channelIndexing.collapse(uu_PhikMean));
    
    scalarField vv_CkMeanValues(channelIndexing.collapse(vv_CkMean));
    scalarField vv_DkMeanValues(channelIndexing.collapse(vv_DkMean));
    scalarField vv_EpikMeanValues(channelIndexing.collapse(vv_EpikMean));
    scalarField vv_TkMeanValues(channelIndexing.collapse(vv_TkMean));
    scalarField vv_PkMeanValues(channelIndexing.collapse(vv_PkMean));
    scalarField vv_PikMeanValues(channelIndexing.collapse(vv_PikMean));
    scalarField vv_PhikMeanValues(channelIndexing.collapse(vv_PhikMean));
    
    scalarField ww_CkMeanValues(channelIndexing.collapse(ww_CkMean));
    scalarField ww_DkMeanValues(channelIndexing.collapse(ww_DkMean));
    scalarField ww_EpikMeanValues(channelIndexing.collapse(ww_EpikMean));
    scalarField ww_TkMeanValues(channelIndexing.collapse(ww_TkMean));
    scalarField ww_PkMeanValues(channelIndexing.collapse(ww_PkMean));
    scalarField ww_PikMeanValues(channelIndexing.collapse(ww_PikMean));
    scalarField ww_PhikMeanValues(channelIndexing.collapse(ww_PhikMean));
    
    scalarField uv_CkMeanValues(channelIndexing.collapse(uv_CkMean));
    scalarField uv_DkMeanValues(channelIndexing.collapse(uv_DkMean));
    scalarField uv_EpikMeanValues(channelIndexing.collapse(uv_EpikMean));
    scalarField uv_TkMeanValues(channelIndexing.collapse(uv_TkMean));
    scalarField uv_PkMeanValues(channelIndexing.collapse(uv_PkMean));
    scalarField uv_PikMeanValues(channelIndexing.collapse(uv_PikMean));
    scalarField uv_PhikMeanValues(channelIndexing.collapse(uv_PhikMean));
    
	//2020.5.21 add scalar budget 
	scalarField s_uu_Dk_1MeanValues(channelIndexing.collapse(s_uu_Dk_1Mean));
	scalarField s_uu_Dk_2MeanValues(channelIndexing.collapse(s_uu_Dk_2Mean));
	scalarField s_uu_Dk_3MeanValues(channelIndexing.collapse(s_uu_Dk_3Mean));
	scalarField s_uu_Dk_4MeanValues(channelIndexing.collapse(s_uu_Dk_4Mean));
	scalarField s_uu_PkMeanValues(channelIndexing.collapse(s_uu_PkMean));
	scalarField s_uu_PikMeanValues(channelIndexing.collapse(s_uu_PikMean));
    scalarField s_uu_EpikMeanValues(channelIndexing.collapse(s_uu_EpikMean));


	scalarField ss_T_1_Dk_1MeanValues(channelIndexing.collapse(ss_T_1_Dk_1Mean));
	scalarField ss_T_1_Dk_2MeanValues(channelIndexing.collapse(ss_T_1_Dk_2Mean));
	scalarField ss_T_1_EpikMeanValues(channelIndexing.collapse(ss_T_1_EpikMean));
    scalarField ss_T_1_PkMeanValues(channelIndexing.collapse(ss_T_1_PkMean));	
	
	
    /*
    scalarField epsilonValues(channelIndexing.collapse(epsilonMean));

    scalarField nuMeanValues(channelIndexing.collapse(nuMean));
    scalarField nuPrimeValues(channelIndexing.collapse(nuPrime));

    scalarField gammaDotMeanValues(channelIndexing.collapse(gammaDotMean));
    scalarField gammaDotPrimeValues(channelIndexing.collapse(gammaDotPrime));
    */

    scalarField urmsValues(sqrt(mag(RxxValues)));
    scalarField vrmsValues(sqrt(mag(RyyValues)));
    scalarField wrmsValues(sqrt(mag(RzzValues)));

    scalarField kValues
    (
        0.5*(sqr(urmsValues) + sqr(vrmsValues) + sqr(wrmsValues))
    );


    const scalarField& y = channelIndexing.y();

    makeGraph(y, UMeanXvalues, "U_x_mean", path, gFormat);
    makeGraph(y, UMeanYvalues, "U_y_mean", path, gFormat);
    makeGraph(y, UMeanZvalues, "U_z_mean", path, gFormat);
    
    makeGraph(y, urmsValues, "uu_mean", path, gFormat);
    makeGraph(y, vrmsValues, "vv_mean", path, gFormat);
    makeGraph(y, wrmsValues, "ww_mean", path, gFormat);
    makeGraph(y, RxyValues, "uv_mean", path, gFormat);
    makeGraph(y, kValues, "k_mean", path, gFormat);

    makeGraph(y, pPrime2MeanValues, "pPrime2Mean", path, gFormat);
    
    makeGraph(y, T_0_001MeanValues, "T_0_001Mean", path, gFormat);
    makeGraph(y, T_0_005MeanValues, "T_0_005Mean", path, gFormat);
    makeGraph(y, T_0_025MeanValues, "T_0_025Mean", path, gFormat);
    makeGraph(y, T_0_05MeanValues, "T_0_05Mean", path, gFormat);
    makeGraph(y, T_0_125MeanValues, "T_0_125Mean", path, gFormat);
    makeGraph(y, T_0_25MeanValues, "T_0_25Mean", path, gFormat);
    makeGraph(y, T_0_71MeanValues, "T_0_71Mean", path, gFormat);
    makeGraph(y, T_1MeanValues, "T_1Mean", path, gFormat);
	
	//add T instane
	makeGraph(y, T_0_001Values, "T_0_001", path, gFormat);
    makeGraph(y, T_0_005Values, "T_0_005", path, gFormat);
    makeGraph(y, T_0_025Values, "T_0_025", path, gFormat);
    makeGraph(y, T_0_05Values, "T_0_05", path, gFormat);
    makeGraph(y, T_0_125Values, "T_0_125", path, gFormat);
    makeGraph(y, T_0_25Values, "T_0_25", path, gFormat);
    makeGraph(y, T_0_71Values, "T_0_71", path, gFormat);
    makeGraph(y, T_1Values, "T_1", path, gFormat);
    
    makeGraph(y, T_0_001Prime2MeanValues, "T_0_001Prime2Mean", path, gFormat);
    makeGraph(y, T_0_005Prime2MeanValues, "T_0_005Prime2Mean", path, gFormat);
    makeGraph(y, T_0_025Prime2MeanValues, "T_0_025Prime2Mean", path, gFormat);
    makeGraph(y, T_0_05Prime2MeanValues, "T_0_05Prime2Mean", path, gFormat);
    makeGraph(y, T_0_125Prime2MeanValues, "T_0_125Prime2Mean", path, gFormat);
    makeGraph(y, T_0_25Prime2MeanValues, "T_0_25Prime2Mean", path, gFormat);
    makeGraph(y, T_0_71Prime2MeanValues, "T_0_71Prime2Mean", path, gFormat);
    makeGraph(y, T_1Prime2MeanValues, "T_1Prime2Mean", path, gFormat);
    
    makeGraph(y, tke_CkMeanValues, "tke_CkMean", path, gFormat);
    makeGraph(y, tke_DkMeanValues, "tke_DkMean", path, gFormat);
    makeGraph(y, tke_EpikMeanValues, "tke_EpikMean", path, gFormat);
    makeGraph(y, tke_TkMeanValues, "tke_TkMean", path, gFormat);
    makeGraph(y, tke_PkMeanValues, "tke_PkMean", path, gFormat);
    makeGraph(y, tke_PikMeanValues, "tke_PikMean", path, gFormat);
    
    makeGraph(y, uu_CkMeanValues, "uu_CkMean", path, gFormat);
    makeGraph(y, uu_DkMeanValues, "uu_DkMean", path, gFormat);
    makeGraph(y, uu_EpikMeanValues, "uu_EpikMean", path, gFormat);
    makeGraph(y, uu_TkMeanValues, "uu_TkMean", path, gFormat);
    makeGraph(y, uu_PkMeanValues, "uu_PkMean", path, gFormat);
    makeGraph(y, uu_PikMeanValues, "uu_PikMean", path, gFormat);
    makeGraph(y, uu_PhikMeanValues, "uu_PhikMean", path, gFormat);
    
    makeGraph(y, vv_CkMeanValues, "vv_CkMean", path, gFormat);
    makeGraph(y, vv_DkMeanValues, "vv_DkMean", path, gFormat);
    makeGraph(y, vv_EpikMeanValues, "vv_EpikMean", path, gFormat);
    makeGraph(y, vv_TkMeanValues, "vv_TkMean", path, gFormat);
    makeGraph(y, vv_PkMeanValues, "vv_PkMean", path, gFormat);
    makeGraph(y, vv_PikMeanValues, "vv_PikMean", path, gFormat);
    makeGraph(y, vv_PhikMeanValues, "vv_PhikMean", path, gFormat);
    
    makeGraph(y, ww_CkMeanValues, "ww_CkMean", path, gFormat);
    makeGraph(y, ww_DkMeanValues, "ww_DkMean", path, gFormat);
    makeGraph(y, ww_EpikMeanValues, "ww_EpikMean", path, gFormat);
    makeGraph(y, ww_TkMeanValues, "ww_TkMean", path, gFormat);
    makeGraph(y, ww_PkMeanValues, "ww_PkMean", path, gFormat);
    makeGraph(y, ww_PikMeanValues, "ww_PikMean", path, gFormat);
    makeGraph(y, ww_PhikMeanValues, "ww_PhikMean", path, gFormat);
    
    makeGraph(y, uv_CkMeanValues, "uv_CkMean", path, gFormat);
    makeGraph(y, uv_DkMeanValues, "uv_DkMean", path, gFormat);
    makeGraph(y, uv_EpikMeanValues, "uv_EpikMean", path, gFormat);
    makeGraph(y, uv_TkMeanValues, "uv_TkMean", path, gFormat);
    makeGraph(y, uv_PkMeanValues, "uv_PkMean", path, gFormat);
    makeGraph(y, uv_PikMeanValues, "uv_PikMean", path, gFormat);
    makeGraph(y, uv_PhikMeanValues, "uv_PhikMean", path, gFormat);

    makeGraph(y, s_uu_Dk_1MeanValues, "s_uu_Dk_1Mean", path, gFormat);
    makeGraph(y, s_uu_Dk_2MeanValues, "s_uu_Dk_2Mean", path, gFormat);
    makeGraph(y, s_uu_Dk_3MeanValues, "s_uu_Dk_3EpikMean", path, gFormat);
    makeGraph(y, s_uu_Dk_4MeanValues, "s_uu_Dk_4Mean", path, gFormat);
    makeGraph(y, s_uu_PkMeanValues, "s_uu_PkMean", path, gFormat);
    makeGraph(y, s_uu_PikMeanValues, "s_uu_PikMean", path, gFormat);
    makeGraph(y, s_uu_EpikMeanValues, "s_uu_EpikMean", path, gFormat);

    makeGraph(y, ss_T_1_Dk_1MeanValues, "ss_T_1_Dk_1Mean", path, gFormat);
    makeGraph(y, ss_T_1_Dk_2MeanValues, "ss_T_1_Dk_2Mean", path, gFormat);
    makeGraph(y, ss_T_1_EpikMeanValues, "ss_T_1_EpikMean", path, gFormat);
    makeGraph(y, ss_T_1_PkMeanValues, "ss_T_1_PkMean", path, gFormat);	
	
	
    /*
    makeGraph(y, epsilonValues, "epsilon", path, gFormat);
    makeGraph(y, nuMeanValues, "nu", path, gFormat);
    makeGraph(y, nuPrimeValues, "nuPrime", path, gFormat);
    makeGraph(y, gammaDotMeanValues, "gammaDot", path, gFormat);
    makeGraph(y, gammaDotPrimeValues, "gammaDotPrime", path, gFormat);
    */
